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Abstract 

In this paper we present calculations of period distribution for old accreting isolated 
neutron stars (INSs). 

At the age about few billions years low velocity INSs come to the stage of accretion. At 
that stage their period evolution is governed by magnetic braking and accreted angular 
momentum. Due to turbulence of the interstellar medium (ISM) accreted momentum 
can both accelerate and decelerate rotation of an INS and spin evolution has chaotic 
character. 

Calculations show that for constant magnetic field INSs have relatively long spin periods, 
> 10 4 -10 5 s, depending on parameters of INSs and ISM density. Due to long periods 
INSs have high spin up/spin down rates, which should fluctuate on a time scale about 
few years. 



1 Introduction 

Spin period is the most precisely determined parameter of a neutron star (NS). Estimates of 
other parameters: masses (for isolated objects), radii, temperatures, magnetic fields etc. are 
always model dependent. Because of that it is very important to have a clear picture of period 
evolution as far as this parameter is usually used to determine other characteristics of NSs. 
Here we try to obtain distribution of spin periods for old accreting INSs (AINSs). 

AINSs are now a subject of interest in astrophysics (see Treves et al.i^). Probably few 
candidates are observed by ROSAT (Motch&). 

In the next section we describe the model we use to obtain period distributions and show 
results for the easiest case of "spin equilibrium". Then in the section 3 we present our main 
results for the "non-equilibrium" case and briefly discuss them in the section 4. Details of 
calculations can be found in Prokhorov et al.&. 



2 "Spin equilibrium" 

Previous attempts to calculate typical periods of AINSs were made by Lipunov & PopovB 
and Konenkov & Popovi\ In these papers the authors do not try to obtain distributions: 
only characteristic periods of AINSs are derived. The authors assume that AINSs are in "spin 
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equilibrium", i.e. all AINSs in these estimates have enough time to reach the stage at which 
magnetic braking is compensated by accretion of angular momentum. We use terms "spin equi- 
librium" and "non-equilibrium" in quotation-marks as far as there is no real equilibrium: period 
can significantly fluctuate. But the situation in general is similar to real period equilibrium in 
close binaries (see Ghosh & LambB, LipunovB). 
We start with the following equations: 

Here I - moment of inertia of a NS, uo = 2ix/p- spin frequency, \x = B R\ — magnetic moment of 

1 /3 

a NS, R co = (GM/w 2 ) - corotation radius, k t - constant of order of unity, and $ - turbulent 
torque, < $ >= 0. J is determined as J =Min (viRcVaRa), where v t = 10 6 cms _1 (Ra/Rt) 1 ^ 3 

- turbulent velocity at R = R G , R t = 2 • 10 20 cm, R G = 2GM/v 2 , v 2 = v 2 + v 2 p , v s - sound 

velocity, v sp - spatial velocity, Ra = (/i 2 /2MV GAfj ^ - Alvfen radius, and va = (GM/Ra) 1 ^ 2 

- Keplerian velocity at the Alfven radius. For the most reasonable parameters J = v t Rc- 

Now we have to introduce a kind of "diffusion coefficient", D, because due to interstellar 
medium (ISM) turbulence we have a kind of diffusion in the space of frequencies. This coefficient 

can be approximatelly determined as D = (1/6) (MJ//) Rg/ v (Lipunov & PopovB). Here 

M = 7r Rq p v - is an accretion rate. 

We can determine an average spin frequency in the following way: 



2 

^turb 



cA-^l)/^/ J Je-vm/D^ V{ \ u i) = i w |3 . (2) 



Finally, p tur b = 2Tr/uJturb- For J = v t ■ Rg we can write it as: 

Pturb = 3.9 • lOV^f M^f 7 V 2 /^ 3/ XL b. (3) 

Here v 7 = ^/(lO^ms- 1 ), R t2 20 = R t /(2 ■ 10 20 cm). 

We note strong dependence of Pturb and D on v. If information on p and p is available the 
problem can be reversed, and one can obtain an estimate of the velocity of a AINS as it was 
done for stellar wind accretion in X-ray pulsars by Lipunov & Popovi*. 

Probability plotted in Fig. 1 was calculated as: 

f{w, v, n) oc (/i 2 /GMLD VV y/D , (4) 

and then normalized. Here V and D are functions of v, /i, n. Note, that peaks in Fig. 1 are 
indeed sharp. It means, that if "spin equilibrium" can be reached, it is possible to use just 
one typical value, p± ur b, for each set of parameters as it was done by Lipunov & PopovS and 
Konenkov & Popovtf . 



3 "Non-equilibrium" calculations 

In this section we calculate probability distribution for the "non-equilibrium" case. These 
distributions are calculated for magnetic moments with Gaussian distribution in logarithmic 
scale with central value lg(//o) = 30.06 and a = 0.32 (see Colpi et al. 2001& for details on 
magnetic evoluyion of NSs). Velocities of AINS are taken from the Maxwellian distribution 
with a mean velocity 200 km s™ 1 . 
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Figure 1: Period distributions in equilibrium for different parameters. In each of the four plots 
we vary one parameter: n, v, M, B. Results are normalized to unity at the maximum. 



We solved numerically the following differential equation: 

df/dt = A/uj 2 d(u A f)/du + D/u 2 d/du (u 2 df/du) , A = fi 2 /GMI. (5) 

After initial parameters of an INS are chosen from the distributions described above we 
check if this NS can reach the accretion stage in 10 10 yrs. To do it we calculate time which 
it spends as Ejector: ts ~ 10 9 n~ 1/ ' 2 f 6 /i3 1 yrs, Vq = v/10 6 cm s _1 . We neglect the stage of 
Propeller, as far as for constant field it is much shorter than the stage of Ejector (Lipunov & 
PopovB). 

In Fig. 2 we present a curve for n — 1 cm -3 and selected INS parameters: fi = 2 • 10 29 G 
cm 3 , v sp = 10 km s _1 . In Fig. 3 we show our final results for Maxwellian velocity distribution 
and log-Gaussian magnetic field distribution for two values of the ISM density. 



4 Discussion 

After an INS come to the stage of accretion it is controlled by two processes (see eq.l): magnetic 
spin-down and turbulent spin- up/spin-down. Initially magnetic spin-down is more significant, 
but at some period, p cr , these two processes become comparable. For longer periods an INS 
will be governed mainly by turbulent forces. One can obtain the following formula for p cr : 
p 2 r = (Ati 2 h 2 )/(GMMJ). An INS reach p cr in At ~ 10 5 -10 7 yrs after onset of accretion: 
At = (J v / ^M)/( / uv / MJ). 
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Figure 2: Period distribution for /i = 2 • 
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Results are normalized to unity at the max- 
imum. 



Figure 3: Period distributions for n — 1 
cm -3 and n — 0.1 cm -3 . Results are nor- 
malized to unity at the maximum of the 
highest curve. 



For field decay picture should be completely different (see for example Konenkov & Popovfl). 
AINSs with decayed field can appear as pulsating sources with periods about 10 s and p about 
10~ 13 s/s. As an INS passes through turbulent cells a value and a sign of p will fluctuate on a 
time scale Rc/v sp ~ 11.8 1^/(10 km s -1 ) yr. 

Acknowledgments 

This work was supported by RFBR (01-02-06265, 00-02-17164, 01-15(02)-99310). 
SP and MP thank M. Colpi, A. Treves and R. Turolla for discussions. 

References 

1 Colpi M., Possenti A., Popov S.B., Pizzolato F.: 2001, in 'Physics of Neutron Star In- 
teriors", Eds. D. Blaschke, N.K. Glendenning, & A. Sedrakian (Springer- Verlag, Berlin), 
fcstro-ph/0012591 ) 



2 Ghosh, P., & Lamb, F.K. 1979, ApJ 232, 256 

3 Konenkov, D.Yu., & Popov, S.B. 1997, PAZh, 23, 569 

4 Lipunov, V.M., 1992, 'Astrophysics of Neutron Stars", Springer- Verlag (Berlin) 

5 Lipunov, V.M., & Popov, S.B. 1995a, AZh, 72, 711 

6 Lipunov, V.M., & Popov, S.B. 1995b, Astron. Astroph. Transactions, 8, 221 

7 Motch, C. 2000, |astro-ph/0008485| 



8 Prokhorov, M.E., Popov, S.B., & Khoperskov, A.V. 2001, |astro-ph/ 01 08503 



9 Treves, A., Turolla, R., Zane, S., & Colpi, M. 2000, PASP 112, 297 



